Fermion condensation around a Coulomb impurity in a Weyl semimetal as a 
manifestation of the Landau zero-charge problem 
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A Coulomb impurity placed in an undoped Weyl semimetal spontaneously surrounds itself with a 
cloud of condensed Weyl fermions. We study this system within the Thomas-Fermi approximation. 
We find that the ground-state of the system is electrically neutral and exhibits an experimentally 
accessible Landau zero-charge effect: the impurity charge is screened out at any finite distance in the 
limit of vanishing impurity size. Specifically, we show how in this limit the Thomas-Fermi equation 
for the electrostatic potential transforms into the Gell-Mann-Low equation for the charge. 
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Over forty years ago Abrikosov and Beneslavskii [1| 
predicted the existence of semimetals having points in the 
Brillouin zone where the valence and conduction bands 
meet with a dispersion law that is linear in the wavevec- 
tor. Such systems, nowadays called Weyl semimetals, 
are likely to be realized in doped silver chalcogenides 
Ag2+sSe and Ag2+$Te [2], pyrochlore iridates Ailr^O-j 
(where A is Yttrium or a lanthanide) Q , and in topolog- 
ical insulator multilayer structures [J]. Weyl semimetals 
make it possible to study in condensed matter laboratory 
settings various effects that were originally proposed in 
the quantum electrodynamics (QED) context. In this 
paper we analyze the screening of a Coulomb impurity 
in a Weyl semimetal; this is related to the vacuum insta- 
bility with respect to electron condensation due to the 
electric field a supercharged nucleus [5j. Our main con- 
clusion is that the character of screening in the Weyl 
semimetal resembles Landau's " zero-charge" situation in 
QED In contrast to QED, the material parameters 

of a semimetal make this effect experimentally accessi- 
ble. We emphasize that the physical mechanism of the 
zero-charge effect described below is different from that 
originally discussed in QED 0, @] ■ Literal analogs of the 
Landau zero-charge effect were predicted to exist in Weyl 
semimetals |l| and in superfluid 3 He — A 8]. 

To demonstrate both the similarities and differences 
between the semimetal and QED settings, let us consider 
an electron excitation characterized by a linear dispersion 
law u = vpk, where Vf is the Fermi velocity and k is 
the magnitude of the wave vector. We ask whether the 
electron (or hole) can be bound by an attractive point 
Coulomb impurity of charge Ze. Assuming the electron 
is localized within a spatial range L, its energy can be 
estimated as 

hv F Ze 2 hv F e 2 

r = —rl 1 - Za )> a = T, — (!) 

where k > 1 is the dielectric constant. Assuming the 
usual order of magnitude estimate vp — 10 8 cm/ sec, the 
"fine-structure" constant a can be estimated as a ~ 1/k. 
Unlike the QED case, both Za < 1 and Za > 1 regimes 



are experimentally accessible. Eq.([T|) seems to imply that 
there is a critical value of the dimensionless impurity 
charge (number of "protons") Z C1 ~ 1/a such as for 
Z < Z C1 the system is not bound at all, L — > 00, while 
for Z > Z C1 the energy |T]) is minimized for infinitely 
localized electron L — > 0. The latter corresponds to the 
atomic collapse in superheavy atoms Q. As in the 
case of QED, the collapse does not occur since in reality 
the impurity attraction at small distances deviates from 
the Coulomb law. 

In QED there exists another (truly) critical charge 
Z C2 170 > Z C1 , for which the total energy of the pro- 
duction of an electron-positron pair vanishes and the vac- 
uum becomes unstable with respect to pair creation; the 
positron repelled by the nucleus escapes to infinity while 
the electron remains near the nucleus [9j . In a semimetal 
the counterpart of Z C2 can be estimated based on the 
observation that due to the zero band gap there is no 
threshold to creation of electron-hole pairs. Indeed, let us 
consider initially bare impurity of dimensionless charge 
Z and ask whether the system energy can be lowered by 
condensing single electron and letting the hole escape to 
infinity. The effect of electron condensation can be un- 
derstood in terms of measurable charge of the system as 
seen at large distances from the impurity. If < Z < 1, 
the measurable charge Z — 1 has a sign opposite to that 
of the bare impurity, and the system energy in this case 
is proportional to (1 — Z) 2 . This would be lower than the 
energy of the bare impurity (proportional to Z 2 ) only if 
Z > Z C2 = 1/2. We conclude that for physically rele- 
vant Z > 1 the condensation of Weyl fermions always 
takes place; interparticle interactions must be accounted 
for when the condensate contains more than one fermion. 

In what follows we will determine the ground-state 
properties of the Weyl semimetal in the presence of the 
Coulomb impurity as a function of Z and a. The phys- 
ical arguments given above suggest that the electrons 
of the "vacuum" (conduction band) spontaneously con- 
dense around the impurity while the holes leave the phys- 
ical picture; the properties of the electron cloud vary with 



2 



Z and a and are determined by the interplay of attrac- 
tion to the impurity (promoting electron condensation), 
and electron-electron repulsion and kinetic energy cost 
of localization of the electrons (which limit the electron 
condensation). For Z> 1 a large number of electrons is 
expected to be present. In this limit we can use Thomas- 
Fermi (TF) theory 10] to compute the properties of the 
system. This approach parallels that employed in QED 
11] to understand the physics of supercharged Z » Z C2 
atoms. Like its more familiar non-relativistic counterpart 
[Toj | , the TF theory developed below becomes increasingly 
accurate asZ^oo and is expected to produce reason- 
able predictions even at Z ~ 1 (l2| . 

The energy of the system can be written as a functional 
of the electron number density n(r), which is assumed to 
be slowly varying with position: 



where 
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The first term gives the kinetic energy of free Fermi-gas 
with linear dispersion relation, and g counts the num- 
ber of Weyl nodes within the first Brillouin zone. For 
example, g = 24 in pyrochlore iridates [3[ and g — 2 in 
a topological insulator multilayer [4]. For simplicity we 
assume an isotropic dispersion characterized by a single 
Fermi velocity vf- The second term is the interaction of 
the electrons with the impurity. The external potential 
<^ea;t(r) is a pseudopotential that represents the pertur- 
bation of the semimetal caused by the impurity; even 
though (fext is not entirely of electrostatic origin we will 
define A(p ext = —Anen ext /K. We assume that the im- 
purity charge density en ex t(r) is spherically-symmetric 
and localized within a region of size a, so that for r > a 
the potential <p e xt(j) reduces to a purely Coulombic form 
fextif) = Zc/ kt corresponding to a net charge Ze within 
the impurity region. The first two terms of the energy 
functional @, if applied to the question of binding of a 
single electron in the presence of point impurity (a = 0), 
lead to an estimate analogous to Eq. ([1} . The last term of 
the energy functional © describes the Coulomb repul- 
sion of the electrons. 

When the number of condensed electrons is large, their 
discreteness can be neglected; then the electron number is 
not constrained and adjusts itself to minimize the energy 
functional Its variational minimization leads to an 
integral equation 
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is a dimensionless parameter that characterizes the rel- 
ative strength of electron-electron interactions and zero- 
point motion represented by the first and third terms in 
Eqs.© and (O, respectively. The weakly- interacting or 
quantum-dominated regime corresponds to the A <C 1 
domain while the classical or strongly interacting limit 
is recovered in the A ^> 1 case. A wide range of values 
for A are experimentally accessible. Taking in Eq. <j3j) the 
r — > oo limit gives (for A ^ 0) a relationship 



n(r)dV = Z 



(5) 



i.e. the semimetal succeeds in giving complete screen- 
ing of the impurity charge; the Coulomb impurity in the 
Weyl semimetal resembles an artificial atom in vacuum 
except that it is made of Weyl fermions. 

When the interactions are weak A <C 1, the screening is 
weak, too. Then neglecting the last term in Eq.Q seems 
legitimate and we find 
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Inside the impurity region where (p ex t — Ze/ na we have 
n ~ XZ 3 /a 3 thus implying that the number of condensed 
electrons residing at r < a, is estimated as XZ 3 . The 
latter is much smaller than the total number Z of con- 
densed electrons provided XZ 2 <C 1 (thus tightening the 
initial A « 1 constraint). Outside the impurity region 
Eq.(j6]) becomes 
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Thus number of condensed electrons confined between 
concentric spheres of radii a and r > a is given by 



N(r >a) = XZ 3 In ■ 



(8) 



Eqs.© and © clearly contradict the normalization con- 
dition ([5]) and the resolution is that there exists a spatial 
scale R beyond which the XZ 2 <C 1 perturbation theory 
and thus Eqs.(J7J) and fail. This scale could be called 
the cloud size; it can be estimated from the condition 
N(R > a) ~ Z with the result R(X) ~ a exp(const/ XZ 2 ) 
(a more precise expression will be given below). 

Eq.j3|), like its non-relativistic counterpart [l0(, can be 
converted into a differential equation by application of 
the Laplacian A: 
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Comparing Eq.© with the Poisson equation for the full 
self-consistent electrostatic potential ip (that combines 
the impurity potential ip ext and electrostatic potential of 
the electron cloud) shows that 



e / Awn(r) \ 



(10) 



The weak-screening XZ 2 -c 1 analysis conducted above 
corresponds to the approximation tp = ip ex t- In terms of 
the physical potential (|TU|) Eq.© becomes 
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For k — 1 and g = 2 (spin-1/2 fermions in vacuum) 
Eq. (ITT1) reduces to the ultra- relativistic limit of the QED 
TF equation [ll|. 

It is convenient to replace ip with the function £ 
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which, via Gauss's theorem, is related to the dimension- 
less measurable impurity charge at position r (defined as 
the total charge inside the sphere of radius r) as: 
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= C{x)-x— = Q(l)-Q(i), i = ln 7 (13) 
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Substituting (1 1 2p into ([9} we obtain the equation 

C"(0 - C'(0 = = -Wa 3 n e;ct + AC 3 (14) 



Substituting (|T2l into (|T0|) we find for the electron den- 
sity the expression of the form 
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Since n must decay faster than 1/r 3 (see Eq.©), the 
solution to Eq. (H"4"[) is subject to the boundary condition 
C(oo) = and to the normalization condition implied by 
Eqs.lJTSj! and ©. 

The weak screening XZ 2 <C 1 analysis can be further 
extended by treating the cubic term of (Q~4|) perturba- 
tively. Then the lowest order solution is £ = Z outside 
the impurity region. The next order gives for r > a 

( = Z(l-XZ 2 l) = Z(l-XZ 2 ln-), AZ 2 Z<1 (16) 

a 

We observe that the expression for the measurable charge 
(|13|l computed with the help of Eq. pl)]) agrees with 
Eq.© to logarithmic accuracy. Then to logarithmic ac- 
curacy the perturbative expression (|16|) may be regarded 
as a charge itself; like Eq.© it is applicable if XZ 2 l -C 1 
i. e. within the Weyl cloud. Eq.(fT6|) tells us that even 
within the cloud, the physical electrostatic potential ip 
and the density of condensed electrons n decrease with 



r slightly faster than Ze/nr and XZ 3 /inr 3 , (see Eq.©), 
respectively. We also observe that Eq. tfTB"]) resembles 
Dirac's QED formula [HI which predicted a logarithmic 
momentum dependence of the physical electron charge 
due to the effects of vacuum polarization. 

As the interactions parameterized by A © increase, 
the role of quantum-mechanical effects weakens. In the 
classical A = oo limit the screening problem has a very 
simple solution: n(r) — n ext (r) and ip = (implied by 
Eqs.© or © and (fTO])). This means that the screening 
response of the semimetal to the presence of the impu- 
rity is that of an ideal classical conductor: the system is 
locally neutral and the electric field is zero everywhere. 

When n ext (r) is a smooth function and A ^> 1, Eqs.© 
and (|10[) imply 



n(r) 



t (r) 



1 



(16n 2 Xy/ 3 



An^ 3 (r) + 



<p(r) 



ATm ext (r) 
X 



1/3 



(17) 



(18) 



We see that for A large but finite the screening is nearly 
complete and a small deviation (of order A -1 / 3 ) of the 
electron number density from n ext (r) is responsible for 
the potential given by Eq. (fT8|) . This imbalanced charge 
is localized in the region of maximal change of the im- 
purity charge distribution n ext {v), i.e. near the impurity 
boundary; the same applies to the residual electric field 
These conclusions resemble earlier results 



11] 



-Vip. 

In practice the strong screening solution ((TT]) and (fi"8|) 
is only relevant inside of the impurity region because the 
impurity boundary has atomic width; a physically ap- 
propriate model of the impurity charge distribution is a 
function n ex t(r) that is strictly zero everywhere outside 
the impurity. For example, for a uniformly charged ball 
model of the impurity Eq.(|17p predicts that n — for 
r > a. However this can be only true in the A = oo 
case; a better approximation requires a solution of the 
non-linear equations © or (fTTj) that does not rely on the 
linearization approximation about n — n ext (r). Eq. (fT8"l) 
allows us to estimate the residual charge within the im- 
purity region to be of the order Ze/ {XZ 2 ) 1 / 3 <C Ze, thus 
implying that most of the condensed electrons are inside 
while only a small number of the order Z / (XZ 2 ) 1 / 3 is 
outside the impurity region. 

This leads to a physical picture of screening as follows: 
*In the weak-screening limit, XZ 2 <C 1, nearly all con- 
densed electrons are outside of the impurity region; of 
the order Z of them reside within the Weyl cloud of size 
R(X) ~ a exp(const/XZ 2 ). Within the cloud, the con- 
densed electron density decays approximately as 1 /r 3 as 



described by Eqs.fTB]) and (fT6j) . while the electrostatic 
potential given by Eqs. (fT2|) and (fTHl) is nearly the same 
as that of the bare impurity, meaning that the screen- 
ing charge is having small effect. Outside of the cloud 
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the electron density and potential decay faster than 1 /r 3 
and 1/r, respectively, in a way yet to be determined. 

*As XZ 2 increases, more electrons invade the impurity 
region and the size of the cloud shrinks. This picture is 
consistent with what we already know about the regime 
of strong screening XZ 2 ^> 1: For XZ 2 = oo the screen- 
ing is complete and all the condensed electrons reside 
within the impurity region; the electric field is zero ev- 
erywhere. For finite XZ 2 S> 1 the screening is nearly 
complete; the majority of the condensed electrons re- 
side within the impurity region, while a small minority is 
found outside. Based on these results we argue that no 
qualitative changes occur in the screening response as the 
system evolves between the regimes of weak and strong 
screening. Below we put this assertion on a solid ground 
by showing that for any A the large-distance screening 
response is universal (in the strong-screening limit this 
will include the entire outside region). 

For I — ln(r/a) > 1 we can neglect in Eq.(fT4|) the 
second-order derivative term C"(0 compared to ('(I) (the 
alternate choice |C"(0I ^ IC'(0I leads to an inconsis- 
tency). Then for large I we have C(0 oc and 
|C(0I > IC(0Ij so that in the / > 1 limit 0) (see 
Eq. ([T3")) ) may be called the enclosed charge. Well out- 
side the impurity region I ^S> 1 Eq. (fT4"|) simplifies to 
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This equation is mathematically identical to the Gell- 
Mann-Low equation 0, [H{ for the physical charge in 
QED [6] reflecting the effects of vacuum polarization. 
Eq. (fT9| exhibits the phenomenon of "zero charge" - no 
matter what the "initial" value of C is, the system "flows" 
to the zero charge fixed point C = as I — > oo, i.e. the 
impurity charge has been completely screened. Since 
/ = ln(r/et), the I —¥ oo limit means either the point 
impurity limit a — > and arbitrary distance r or fixed 
impurity size a and r — > oo. Eq.(fT9"|) can be integrated 
with the result 
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(20) 



where Co is an integration constant determined by the net 
charge within the impurity region. At sufficiently large 
length scales, / — ¥ oo, the initial value Co is unimportant 
and C is independent of Z, for any value of A. 

Substituting Eq. ([2U)) into the expressions for the po- 
tential (fT2")l and electron density (|15l) we find our main 
results 
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applicable in the I = ln(r/a) ^> 1 limit. These equa- 
tions can be viewed as interpolation formulas for arbi- 
trary r a, A, and Z, with Co having the meaning of 
the net dimensionless charge within the impurity region. 
The logarithmic terms in the denominators of Eqs. (|21[) 
and (|22p are relevant at the scales r exceeding 



R(X) ~ ae 1 ^ 



(23) 



which is the size of the Weyl cloud. 

The total number of the condensed electrons confined 
between the spheres of radii a and r > a computed with 
the help of Eq. ([22|) 
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(24) 



is finite and equal to Co as r — > oo. 

In the weak screening limit XZ 2 <C 1 when 
Co = Z, the size of the Weyl cloud (|2"3")) becomes 
R(X) ~ aexp(l/2A2' 2 ), sharpening the previous esti- 
mate. Moreover Eqs. (|20| and (|24[) include the pertur- 
bative XZ 2 hi(r/a) <C 1 results, Eqs. ([TBI) and ([HJ, respec- 
tively. We also observe that in the weak screening limit 
the Weyl cloud described by Eqs. (|2"Tl) and (|2"2"|) is an ex- 
tremely weakly localized object: only about 30% of the 
condensed electrons are found within the range (|23|) . 



In the strong screening limit XZ ^> 1 when Co — 
Z/iXZ 2 ) 1 / 3 < Z, the size of the Weyl cloud (J23J) ap- 
proaches the impurity size. This means that asymptotic 
universal regime described by Eqs. (|2"Tj) and (j2"2"j) is the 
solution to the strong-screening problem practically ev- 
erywhere outside the impurity region. 
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